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SUMMARY 


Values of coefficients for defining the effectiveness of integral 
stiffeners in resisting shear deformations of the plate of which they 
are an integral pant are presented for a wide range of proportions of 
rectangular stiffeners with circular fillets. The coefficients are 
evaluated by the use of an electrical analog computer. Formulas are 
given in which these coefficients may be employed to calculate the 
elastic constants associated with the twisting and shearing of integrally 
stiffened plates, either directly, as in the case of simple longitudinal 
or transverse stiffening, or throijgh the intermediate evaluation of pre- 
viously defined shearing-effectiveness parameters, as in the case of more 
complicated stiffener patterns . The fillet radius is shown to contribute 
appreciably to the degree of penetration of the stresses from the skin 
into the stiffeners. Thus, through the use of suitable combinations of 
rib proportions and fillet radii, simple longitudinal or transverse 
integral stiffening can be made to contribute to the overall shear 
stiffness of the plate-stiffener combination. 


INTRODUCTION 


The effectiveness of integrally stiffened plates, as demonstrated 
in references 1 to h, is in part due to the fact that the stresses in 
the skin of the plate are conducted into the integral stiffeners. Even 
under simple loadings, a complicated stress distribution within the 
cross section is produced, and in consequence the evaluation of the elas- 
tic constants for the plate- stiffener combination is difficult. 

Form\ilas for the thirteen elastic constants of plates with integral 
stiffeners and a method for analytically obtaining upper and lower limits 
on the parameters of the formulas are presented in reference 4. In some 
cases the differences between the constants calculated by the upper- 
and lower -limit assumptions are substantial. For example, upper -limit 
shear stiffnesses calculated in reference 4 were from 12 percent to 
52 percent greater than the calculated lower-limit stiffnesses. 
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In the present paper a more refined analysis is made of the shearing 
effectiveness of integral stiffening. The approach used is that of 
imposing a quasi-shear deformation upon a repeating element of a plate 
vith simple longitudinal or trainsverse stiffening (see fig. 1) and then 
solving the eqxiations of elasticity associated with the imposed deformation. 
The method of solution is similar to the method used in solving the torsion 
problem of pages 258 to 265 of reference 5. This approach requires the 
solution of Laplace's equation over the cross section of the repeating 
element as shown in figure 2. Solutions for a wide range of proportions 
of rectangular stiffeners with circiolar fillets were obtained with a 
General Electric Analog Field Plotter (ref. 6) which was modified by 
the MCA to suit the needs of this particular problem. This modified 
field plotter is similar In operation and principle to an electrical 
analog computer described in reference 7 • 

Results of this analysis are presented in the form of tabulations 
and curves giving coefficients from which the shearing effectiveness of 
the integral stiffener may be evaluated. These coefficients may be 
used with the formulas of reference k for the calculation of the plate 
elastic constants; if the shearing stiffness of a plate with simple 
longitudinal or transverse stiffening is required, however, it may be 
determined more directly from the given coefficients throu^ the use of 
formulas presented herein. 


SYMBOLS 

Plate Dimensions 


A^ area of perpendicular cross section of rib, sq in. 

bg length of repeating element of integrally stiffened plate, in. 

by height of rib above plate, in. 

H total height of rib and plate, tg + by, in. 

h z-distance from y-axis to boundary cdef of specimen (fig. 2) 

Ig length of plate between fillet and end of repeating element, 

bs - 

2 

Zy height of rib above fillet radius, by - ry, in. 


NACA TN 3443 


3 


■^S 


D 


xy 


xy 


K 

K' 

^xy 



radius of fillet, in. 
thickness of skin or plate, in. 
thickness of rib or web, in. 


Forces and Elastic Constants 

coupling elastic constant associated with coupling between 
twist and shear and defined by equations (l) and (2), in. 

twisting stiffness relative to x- and y-directions defined 
by equation (l), in-lb 

twisting stiffness relative to x- and y-directions defined 
by equation (3), in -lb 

shear modulus of material, psi 

shear stiffness of plate in xy -plane defined by equation (2), 
Ib/in. 

shear stiffness of plate in xy-plane defined by equation (4), 
Ib/in. 

torsion constant defined by equation (27) 

torsion constant for sections shown in figure 6(a) 

intensity of resultsint twisting torque, lb 

Intensity of resultant shearing force acting in plane z = t^2, 
Ib/in. 

coupling elastic constant associated with coupling between twist 
and shear defined by equations (3) and (4), lb 


Special Symbols Used in Shearing-Effectiveness Analysis 
coupling coefficient defined by equation (A30) 
coupling coefficient defined by equation (A25) 
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coupling coefficient defined by equation (A26) 


xy 


■yx 


J 

P 


■"xy 

a’ 

3 ’ 

3 


intensity of resultant shearing force acting in y-direction 
in plane z = ts/2, Ib/in. 

intensity of resultant shearing force acting in x-direction 
in plane z = tg^P, Ib/in. 

coefficient of twisting stiffness defined by equation (6) 

boimdary value of stress function or boundary value of 
electrical potential field (fig. 2) 

coefficient of shearing stiffness defined by equation (A22) 
and determined by equation (A24) 

integral of stress function defined by equation (A55) 

integral of stress fxmctlon defined by equation (A 56 ) 

intensity of resiiltant twisting torque acting on planes per- 
pendicialar to x-axis, lb 

intensity of resultant twisting torque acting on planes 
perpendicular to y-axis, lb 

coefficient used in reference 4 to locate effective centroid 
of part of rib for resisting twisting deformation 

coefficient vised in reference 4 to define effectiveness of rib 
in resisting shear deformation 

coefficient used in reference 4 to define effectiveness of rib 
in resisting stretching in its -transverse direction 

magnitude of pure shear distortion Imposed upon repeating 
element of integrally stiffened plate (fig. 3), In. 


General Symbols 

u,v,w displacements in x-, y-, and z-directions, respectively, in. 

orthogonal coordinates j z measvired normal to plane of plate, 
and X and y measured in plane of plate 


x,y,z 
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^xz^^yz 


^xy 


O'x^ Oy> ^z 

^xy’^xz^''’yz 




components of shear strain 

average r^_ over length of repeating element defined 
•xy 

by equation (A21) 
components of normal strain 

con 5 )onents of normal stress, psi 
componen-ts of shear stress, psl 

transformed orthogonal coordinate system defined in 
equations (22 ) 

stress function defined by equations (a 6) 
warping function defined by equation (A1) 


STATEMENT OF PROBLEM 


The force-distortion relationships for the twisting and shearing of 
rectangular orthotropic integrally stiffened plates having their axes 
of principal stiffness parallel or perpendlciilar to the sides of the 
plate (figs. 1 and 2) may be written, as noted in reference h, in two 
forms as follows; 


or 


;v2 

o w 




S^w , ^xy 
^xy k Sx 3y Gk 


Sx Sy Djjy ^ 


^xy 

y = 2TM + y; 

^xy xy 


( 1 ) 

( 2 ) 


( 3 ) 

i^) 


6 


NACA TN 


in which the elastic constants D , G , and T are related to D 

xy' xy' 

Gj^, and Cj^ by the following formulas: 

= 2Di. 




xy 


xy 


T = 




®k + ^k ^k 


2D,. 


(5) 


These elastic constants can be evaluated in terms of the three coef- 
ficients j, q, and c which express the effectiveness of the stiffeners 
in resistance to twisting, to shearing, and to coupling between twisting 
and shearing, respectively. For plates with simple longitudinal or 
transverse integral stiffeners, the equations for the elastic constants 
in terms of these coefficients are determined as follows: 


I>k = 1 

(6) 

Gk = 

(7) 

^k = - ''■^3 

(8) 

“xy • GJts’ 

(9) 

G jqtq 

G = ^ 

(10) 

c 

.. 2 

(11) 


The determination of the coefficient j is presented subsequently in the 
text and that of the coefficients q and c, in the appendix. 
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Expressions for these elastic constants have also been derived in 
reference 4 for plates with stiffeners in a variety of patterns, where the 
effectiveness of the stiffeners for resisting twisting and shearing is 
expressed in terms of two parameters a' and p'. A method for the 
evalixation of a' was given in reference 4, based on the work of refer- 
ence 8, but no basis for the evaluation of p' was then available. 

(The suggestion was made in reference 4 that p' be assumed equal to p. 
Subsequent experimentation has shown that for a wide range of proportions, 
p = 7/8 p'.) Values of a' and p' may now be determined, however, 
in terms of J and qj thus 


where 









( 12 ) 

(13) 


(14) 


Equations (l2) and (l5) are derived in the appendix. 

The problem considered in the present paper is the evaluation of j, 
q, and c . Actual values of these coefficients are obtained herein for 
only rectangular stiffeners with circular fillets, but the methods of 
analysis are applicable to stiffeners of any cross section. 


PROCEDURES FOR DETERMINATION OF COEFFICIENTS 
Determination of j 


References 8 and 9 give torsion constants for sections such as those 
illustrated in figures 1 and 2. The coefficient j used in the evalu- 
ation of the elastic constants of integrally stiffened plates is related 
to the conventional torsion constant K as follows: 


K 



( 15 ) 
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Check tests of the twisting stiffness of a few plates of different 
proportions have been made which confirm the results of references 8 
and 9- No further analysis of the twisting stiffness of integrally 
stiffened plates was therefore made. Further Information on the evalu- 
ation of j is given in the section entitled "Results and Applications." 


Determination of q. and c 

A quasi -shear deformation was imposed on the repeating element of 
the integrally stiffened plate as shown in flgme 3 , and the stress 
resultants required to produce this deformation (fig. 4 ) were found. 

The details of this analysis are presented in the appendix and the 
resulting equations are found to be 

Nxy = <lGtg7xy ( 16 ) 

and 

“xy - - (IT) 


where q and c are given by 


la 

2 :^+ Si 
tg 


(18) 


^ 1 ^ Ji! ^ ^ 2(1 . 

tg 2 tg tg tg V 4 / \tgy 

c = : 


S2 - I Si 


l*sAs 


( 19 ) 


in which bg, by, Iq, tg, and ty are dimensions of the plate 

(see fig. 2), and Si and are integrals of a stress function 0 

defined in the appendix. The integrals are 


S 


1 



0 


z=t 


dy 

s 


2J 
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and 



h 


S 


0 dz dy 


In order to evaluate these integrals, use was made of an electrical 
analogy. As is shown in the appendix, the stress function 0 must 
satisfy Laplace’s equation. 


The application of Ohm's law to a thin conducting sheet of material 
(see ref. j) shows that a function V, which describes an electrical 
potential field in the sheet, must satisfy Laplace's equation. 


Accordingly, a potential field was set up in a conducting sheet over a 
shape related linearly to the shape of the cross section under consider- 
ation (see fig. 2) by electrically duplicating the boxmdary conditions 
on the stress function. A self-balancing potentiometer was used to 
measure values of the potential over the conducting sheet as shown in 
figure 5* The desired integrals and Sg were then computed by 

numerically integrating the potential readings over the cross section. 

The potentiometer used was a modified General Electric Analog Field 
Plotter (ref. 6), and the conducting sheet was Type L Teledeltos paper. 
This paper was sufficiently homogeneous to give good results with the 
large-scale cross sections used, but an adjustment was required to take 
into account a 7- percent deviation from the maximum resistance which 
existed between the directions of principal resistance. This directional 
property was corrected for by distorting the figure according to the 
transformation 



( 20 ) 





( 21 ) 


z = ^ 


( 22 ) 
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where Ry and are electrical resistances in the y- and z-directions, 

respectively. The actual potential field in the sheet is described by 


^2 p ^2 

5 V ^ ^ 


(23) 


Thus, a potential value measvured from the distorted figure is equivalent 
to a potential at a corresponding point in an undistorted figure on 
uniform, nondirectional conducting paper. 

As a check on the overall accuracy of the procedure, values of the 
integrals were calculated for a typical cross section by the iterative 
procedure of reference 10 and compared with those measured with the 
electrical analog. The difference between the two values was less than 
1 percent. 


RESULTS Airo APPLICATIONS 


The values of j, q, and c may be used directly in equations (6) 
to (ll) for the calculation of the elastic constants associated with the 
twisting and shearing of plates with simple longitudinal or transverse 
integral stiffening. For plates having combined longitudinal and trans- 
verse or symmetrically skewed ribbing, values of j and q may be used 
in equations (l2) and (13) to calculate corresponding value of the param- 
eters a' and p' for use in the elastic-constant formulas of refer- 
ence 4. 


Evaluation of Coefficients 


Evaluation of j . - Check tests performed in the Langley structures 

research laboratory in conjunction with this investigation have shown 

that the method of reference 8 gives accinrate results only when 

^w/^S ^ following formulas, which can be derived 

by use of reference 8, should therefore be used only when ry^tg % 1 
and ty/tg ^ 2: 



1 i 

2 bg 3 





( 24 ) 
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where 



(25) 


and 


_d 

tc 


S 



(26) 


The results of reference 9 were foimd to agree closely with the 
check tests previously mentioned. Reference 9 presents torsional 
constants (designated herein as K') for the configuration shown in 
figiu-e 6(a), in which = "tw- results of refer- 

ences 8 and 9 and extrapolations of each of these sources are presented 
graphically in figure 6(b) for integrally stiffened plates of the pro- 
portions shown in figure 6(a). The extrapolations were partly guided 
by the experimental data from the check tests. 

Inasmuch as the curves of figure 6(b) apply to only the config- 
uration shown in figure 6(a), the torsional stiffness of any additional 
skin or rib height must be accounted for separately. Since the addi- 
tional skin or rib is remote from the juncture, the stress distribution 
at the junct\ire will not significantly affect the distribution in the 
remote portions of the section; therefore, the additional torsional 
stiffness is approximately that of a rectangular section. The torsion 
constant K for any section is then the sum of K' (the torsion constant 
of the section shown in fig. 6(a)) and the torsion constants of the rec- 
tangular sections; thus 



The relation between K and j is given in equation (15). When exper- 
imental data are available, J may be computed directly as 



Gtg 8 w 
8x 8y 


(28) 
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Evalmtlon of q and c . - The values of the shearing and coupling 
effectiveness coefficients q. and c determined by the electrical analogy 
are given in table I for a wide range of proportions of rectangular ribs 
with circxilar fillets. Also included in table I are the corresponding 
values of the parameter p' of reference 4. 


Values of q 


^3/~^S 

<1 



and 


and c are presented in figures 7 8 plots of 

t)g ry 

— c against — . The ciirves apply only when the 


values of both 




and 




axe equal to or greater than the 


largest values of these two parameters which appear in the table for the 
given values of ry/tg and ty/tg under consideration. These plots 


may be made because, beyond certain limits 


(the maximum values of 




and by^tg appearing in table I for the ry^tg and tyy'tg under con- 
sideration), additional rib height or additional plate length between 
ribs will not affect the stress distribution at the juncture of the skin 
and rib. 


Interpolation may be made by cross-plotting when a set of dimension 
ratios fall within the range of values presented in the table. 


When bg^tg is found within the table but by^tg lies beyond the 

range of the table, the values of q and c appearing under the largest 
value of by/tg for the bg/tg being considered may be used. However, 

the value of 3' must then be computed from equations (13) and (l4) by 
using the actual value of ^y/'tg* When by/tg is fo’und within the table 

but bg/tg lies beyond the range of the table, the value of q for the 
largest value of ^g/fg for "the value of t)y^tg being considered is used 

in equation (l8) to obtain a value of The value of bg/tg for which 

the value of q was foiand must be used in equation (l8) when is 

computed. The value of thus obtained is then used in equation (l8) 

with the required values of tjg/'tg and ^g/tg to obtain the desired 

values of q. This value of q may then be used in equation (l5) to 
compute 3'- A similar scheme is lised to find the value of c. That 
the value of c found in the table corresponding to the value of the by/tg 
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under consideration and the highest value of l>s/tg appearing for that 

value of by/tg, together with the value of S 2 _ previously obtained, 

are used in equation (19) to compute S 2 . The values of 82 ^ and 

the dimension ratios being considered are then used in equation ( 19 ) to 
compute c . 

The ranges of the dimension ratios covered by the table and curves 


are : 

Ratio Range 

Rib thickness/Skin thickness, ly/tg 0 to i)- 

Fillet radlus/Skin thickness, ry^tg 0 to I 6 

Length of repeating elements/Skin thickness, Lg^tg 0 to 00 

Rib height/Skin thickness, b^/tg 0 to «> 


An additional result incidental to the evaluation of q and c 
is that, for values of r^/tg > 1, no shear stress concentration exists 

in the cross section. (Stress-concentration factors were based on an 
average shear stress in the skin at a remote unaffected distance from 
the juncture of the skin and ribs.) No investigation of the stress 
concentration when ^ ^ pursued. 


Illustrative Examples 

In order to illustrate the method of obtaining the effectiveness 
coefficients and their significance, the effectiveness coefficients j, 
q, and c are calculated for the rib proportions used in the tests of 
reference 3 and one variation of that shape. 

The dimensions of the rib cross section used in reference 3 are: 
bg = 1.00 in. tg = 0.05 in. 

ty = 0.10 in. by = 0.20 in. 


The dimension ratios are then 


%/‘ts = 2 

bg/tg = 20 


ry = 0.10 in. 
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ry/tg = 2 


These dimension ratios are considered in the first example. The second 
example considers the same dimension ratios except that the value of 
ry/ts is changed from 2 to 0. 


Example 1.- From figure 6(b), for — 

tg 


2 and 


fw 


2 , 



Therefore, 


K' 


23 . 4t 


4 

S 


From equation ( 27 ) the torsion constant K is determined as 


K = 23.4tg^ + |(l4 - 2)tg^ + j(2 - 2)(2)\g^ = 27.4tg^ 

The effectiveness coefficient j can now be computed from equation ( 15 ) 
as 


li 

27.4tg 

2(20)tg^ 


0.685 


The value of bg/tg is larger than the largest value of ^ 3 /tg which 
appears under ty/tg = 2 and ry/tg =2 in the table, and by/tg is 
equal to the largest value of by/tg which appears under ty/tg = 2 and 
ry/tg = 2. Figures 7 anb 8 are therefore used to obtain values for q 

bg/tg Ig 

and c. From figure 7 , the value of 2 — for i*y/tg = 2 and 

^ S 

ty/tg =2 is found to be 4.0. Hence 

^ - 14 = 4 
<1 
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or 


q = 1.111 


From figure 8, 



is found to be 2.10, or 


c = 0.105 

Using the above values of J, q, and c in equations (6), (7), and (8) 

gives the values of Dj,, G^, and as follows: 

= 42.8 X 10"^G 
= 0.0555G 

Cj^ = - 0.00525 

The value of a' is found from equation (12) by substitution of the 

values of j and q previously found: 


a 




0.3225 


From equation (l4). 




bgtg 20 


4(2) + 2(1 - |)(2)' 


= 0.486 


The value of 3' is then found from equation (13) as 

1.111 - 1 


3' = 


0.486 


= 0.2285 


Example 2 .- In this example, a configuration having the following 
proportions is considered: 

%As “ ^ ^sAs = 

^wAs " ^ ^wAs ° 


16 


NACA TN 3 hh 3 


The value of j may be found by means of equation ( 24 ) or from fig- 
ure 6 (b) and equation (15) to be 

j = 0.456 

From figures 7 8 ^ 0 . c are found to be 

q = 1.031 


and 


c = 0.032 

The values of and are then found from equations ( 6 ), (j), 

and ( 8 ) to be 

Dk = 27.30 X 10 "^G 

&k = 0.0515 50 
Ck = - 0.00160 

From equation ( 12 ), 

From equation (l 4 ). 


The value of p' is then fo\md from equation (l 5 ) to be 

P' = 0.0775 

The values of a' and p' found in these two examples could have 
been used in the formulas of reference 4 to obtain the elastic con- 
stants Djj., G^, and C^. The values of Dj^. and G^^ so obtained wo\old 

be exactly those values obtained by using the effectiveness coefficients j 
and q in equations ( 6 ) and (j). The values of obtained from the 

fonuulas of reference 4 , however, are found to be somewhat greater than 
those computed directly by use of the coefficient c. 


a' = 0.424 




'=s*s 


• OAO 
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This discrepancy, which disappears when ribs of small twisting 
stiffness are considered, arises as a result of the assumption of refer- 
ence 4 that the shearing effectiveness of the ribs can be represented 
by a substitute sheet of zero twisting stiffness. There is no corre- 
sponding discrepancy in the calciilation of Dj^, however, since in refer- 
ence 4 the location of the substitute sheet (measured by a'H) is chosen 
to give the correct value of twisting stiffness for the stiffened plate 
as a whole. If a value of a' is desired which will give the correct 
value of the coupling term that value of a' may be obtained by 

equating the expression of reference 4 for to that of the present 

paper pmrt solving for a' . That procedxare leads to the following expres- 
sion for a' : 



This value of a' would give correct values of C-^ but somewhat con- 
servative values of Dj^. 

Discussion of illustrative examples .- The most significant implication 
of the results of this evaluation of the effectiveness coefficients is 
that relatively small changes in detailed proportions can appreciably 
affect the overall effectiveness of integrally stiffened plates. As dem- 
onstrated in the examples, a change in fillet radius from ~ ^ 

to r^/t 3 = 0 decreased 3' ( 3 ' is a parameter which shows the effi- 

ciency of the rib in resisting shear) by a factor of 0.339 and decreased 
the twisting stiffness by a factor of 0 . 687 - 

More conqslete analysis will be required to evaluate fully the merits 
of large fillet radii and the effect of changes in the other proportions. 
Results of buckling tests of plates with fairly large fillet radii, such 
as those of reference 3, should not be considered representative of the 
results to be expected for similar plates with small fillet radii. 

Suitably proportioned longitudinal or transverse integral stiffening 
can evidently contribute to the shear stiffness of plates; thus integral 
stiffeners may be utilized to contribute to the torsional stiffness as well 
as to the flexural stiffness of wing panels. 
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CONCLUDING REMAEKS 


The evalijation of the shearing effectiveness of integral stiffening 
for a wide range of proportions of rectangular stiffeners with circular 
fillets has indicated that the degree of penetration of stresses from 
the skin into the stiffeners is in part dependent upon the fillet radius. 
Also, for fillet radii greater than the skin thickness, the shear-stress- 
concentration factor has been found to be equal to \mity. Determination 
of the overall structural importance of the fillet radius and the effect 
of changing other proportions require and, on the basis of the large 
changes in stiffness associated with small changes in configuration shown 
in the present study, deserve further investigation. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 1955* 
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APPENDIX 


ANALYSIS OF SHEARING EFFECTIVENESS OF LONGITUDINAL 


OR TRANSVERSE INTEGRAL STIFFENING 


Figlire 5 shows the repeating element of a plate with integral 
unidirectional stiffeners. A quasi-shear state of deformation has been 
imposed upon the element so that the edge af has undergone a pure 
shear translation with respect to the edge be. The problem is to 
determine the stresses necessary to maintain the imposed deformation. 
The ratio of the resultants of the stresses to the magnitude of the 
assumed distortion provides an index to the effectiveness of integrally 
stiffened plates in resisting shearing forces. 


Derivation of Differential Equation Governing Stress Function 

The semi-inverse method of Saint-Venant, as found in reference 5 
(pp. 259-265), is the approach used for this problem. Plane sections 
parallel to the yz-plane before distortion are assumed to have their 
shape preserved, but these planes may warp in the x-direction. This 
warping is the same for all cross sections along the x-axis. Displace- 
ments u in the x-direction of points in cross sections parallel to the 
yz-plane can be defined by a warping function. 


Since the shape of the cross section is preserved, the displacements in 
the y- and z-directions (v and w) , respectively, are 


The components of strain are therefore calculated from the relations 
between strains and displacements as 


u = i(y,z) 


(Al) 


V = w = 0 


(A2) 


Sx - “ 7yz “ 



(A5) 
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The corresponding stresses can then he calculated as 


cTx - - o'z - ‘’■yz 


= 0 


= G ^ 
z 


Txz - G ^ 


T - r 


(AW 


Consider now a stress function 0 = 0(y^z) from which the shear 
stresses Tjjz T'xy obtainable. The equations of equilibrium 

given on page 229 of reference 5 must be satisfied. Only one of the 
three equations of equilibrium is of significance, namely, 

^ ^ = 0 (A5) 

Sz hy 


The stresses t„_ and may then be expressed in terms of the 

XZ Xjr 

stress function 0 = 0(y,z); thus 


' xz 


= M 


'*’xy “ 


hz 


(A6) 


Equating the stresses determined in equations (Ak) to those determined 
in eqioations (a 6) , so as to determine the stresses from a consideration 
of displacements and thereby automatically satisfy con^iatibility of 
strains, yields 


^z Sy 

G ^ 

by hz ^ 


(A?) 


The warping function if may be eliminated from equations (A7) by 
differentiating both sides of the first of equations (A7) with respect 
to y and both sides of the second with respect to z and subtracting 
the second from the first. Elimination of the warping function shows 
that the stress function 0 miist satisfy laplace's differential 
equation 




(A8) 
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g-^j* 0 gg 0 Q ob'b&in©! f*!Toni "tli© solu.'bioii of* 'fcliis d.if*f*©r©n'fci&l ©Q^ufl'bioii s&tisfy 
the conditions of equilibrium and compatibility. 


Determination of Boimdary Conditions 

The stresses normal to the boundaries ab and cdef of the element 
must be zero because these boundaries represent the stress— free surfaces 
of the repeating element. (See fig. 5-) The values of the stress 
function along these two boundaries must therefore be constant in order 
that the stresses normal to these boundaries may be zero. Boun dapr ab 
is arbitrarily set at the constant value of zero, and boundary cdef 
is arbitrarily set at the constant value of P. The physical signifi- 
cance of this choice of boundary conditions may be seen by considering 
the integrally stiffened plate to be a flat plate. Then, 



is the solution to eqimtlon (a 8 ) . Thus, 



(A9) 


(AlO) 


or 

= -^xyts = -I* 

Therefore P is the magnitude of the applied shear force per unit 
length. 

Along the boundaries af and the displacements u are assumed 

to be constant. The shear strain 7^z those boundaries is given 

y = ^ + ^ = 0 (A12) 

The shear stress along those boundaries is therefore calculated 

to be zero, or 

M = T = 0 (A15) 

8y 

along boundaries af and The assun 5 )tion that = 0 along 

boundaries af and be could have been made from consideration of 
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the fact that the houndaries are lines of symmetry for the stress 
function, and therefore = 0 at these lines of syrametry. 

The problem has now been reduced mathematically to solving equa- 
tion (a 8) subject to the boundary conditions: 


0=0 along ab 


0 = P along cdef 
^ = 0 along af and be 


(A14) 


which are included in figure 2. 


Determination of Resultant Forces 

The resultant forces (see fig. 4) necessary to maintain the assumed 
state of deformation can now be determined in terms of the stress 
function 0 = 0(y,z). When these forces and the distortions produced 
by them are determined, the elastic constants can be obtained in terms 
of the stress function. 


In the plane y = bg/2, the shear stresses may be resolved to 
a resToltant shearing force Fy^ per unit length, acting at the midplane 
of the skin, and a torque Ty^^ per unit length. Equations (a 6) and 
the boimdary conditions (eqs. (Al4)) are used to determine the forces 
in the plane y = bg/2 as follows: 
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In a like manner, the resultant forces on the yz-plane are obtained as 

-'h 


1 


.bg/2 “0 


dy = -P 


(ait) 


and 


_ ^ r'-s/s r- 

^ ‘ bS J.bs/2 ^ 


'’’xyy^ ^xz^ 


dz dy 


= Ptsf- + 2 


2 r^s/2 rn 


1 


2 bsts Phsts ^_-bg/2 ^0 


0 dz dy + lV=l>s/2 


dz 


(A18) 


where h is the z-distance from the y-axis to the boundary cdef. 


Derivation of Formula for q 

The relative shear displacement 5 of plane af at y = bs/2 
with respect to plane be at y = -bs/2 at any value of z between 
0 and ts is given by 

pbg/2 

S = / Txv'^ (0 a z = "ts) (AI 9 ) 

^-bs/2 


From this equation, eqioations (a 6) , and the stress-strain relation- 
ship Txy = G7xy, S becomes 



(0 < z < ts) (A20) 


The average shear strain over the length of the repeating element 
is given by 



1 

Gbgts 



$^z=tg^y 


(A21) 
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From equation (AI 7 ) , = Fy^ = -Pj therefore, each of these 

two shear forces can he replaced hy the more conventional notation N^y. 
The coefficient q, which represents the effectiveness of the integral 
stiffener in resisting shear, can now he defined as follows: 


1 = 


Gts/xy 


(A22) 


Substitution from equation (AI 5 ) or (AI 7 ) for N^^y and from equa- 
tion (A21) for 7 xy in equation (A22) yields the coefficient q in 
the following form: 


q = 



(A25) 


Noting that 0 = P when z = ts and ry + ^<|y|'$ T'S + ^ 

(see fig. 2 ) and expressing the dimensions as dimensionless ratios 
permit q to he obtained in the form 


q = 



1 

Pts 


~bs/ts 


tw 



0z=ts<iy 


(A24) 


Derivation of Formula for c 

The resultant torques Tjjy and Ty^^ necessary to maintain the 
assumed state of deformation are defined in terms of two coeffi- 
cients crp and Cm as 
•^xy -^yx 

Txy = -CT^t3Njjy 
■^Tyx'^S^xy 

The torques are negative if Nj^y is positive. This can he seen hy 
comparing Fyx of equation (AI 5 ) and Ty^ of equation (AI 6 ) . From 
equations (AI 7 ) , (AI 8 ) , and (A25), cm is obtained as 


(A25) 

(A26) 
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— - — / / 0 dz dy 

^ - 1 ) 3/2 '^0 


,bg/2 oh 



,dz (A27) 



1 

Cm = — - 

■^yx 2 



i^y=hs/2' 


idz 


(A28) 


In order that these results may he incorporated in flat-plate theory, 
the twisting moments on adjacent sides of a repeating element must be 
equivalent. The distributed moments Tjjy and may be replaced by 

concentrated lateral forces T^y and Ty.^ at the corners of the plate 
as is done in reference 11. The resultant torques and Tyjj. are 

then replaced by where 


This system is statically equivalent to the actual system. It therefore 
produces essentially the same distortions as the act;ial system except in 
regions at the edges of the plate comparable in width to the thickness 
of the plate. This result follows from Saint-Venant 's principle (see 
ref. 5, p. 55). The resultant torque is therefore defined in 

terms of a coefficient c as 




(A29) 


- -ctgWj^ 


(A50) 


where 




From equation (A20) the following relationship exists: 



(0 < Z < tg) (A52) 
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Rearranging equation (A52) and integrating over the thickness of the 
skin yields 



tg 

2 '^-bs/2 


0.=tg<iy 


(A33) 


Tlius, in finding the value of a double integral of the 0-surface over 
the area of the element^ only the integral over the attached stiffener 
need he founds since the integral over the flat-sheet part of the 
element is found from just the line integral along the line z = tg 
over the length of the element. 


Substituting equation (A55) into equation (A5I ) y noting that 0 = P 
when z = tg and 1^1 5 nondimens ionali zing 
yields 


c 



1 ^ 
2 ts 


+ 


^ % 
ts ts 




0 dz dy - 


1 

2Ptg 



0z=ts^y 




(A54) 


Thus, the forces required to maintain the assumed state of deforma- 
tion are known in terms of two integrals : 



(A55) 


and 



(A56) 


These two integrals are evaluated for the cases considered herein by 
means of an electrical analog corr5)uter. 


Determination of Elastic Constants Gry^y and in 

Terms of q, and c 

The elastic constant is defined in terms of the coefficient j 

in equation (6) where j is obtained by the methods discussed in the 
main text. 
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The state of deformation assumed in the determination of the 
coefficients q and c is characterized by 

= 0 (A57) 

Sx Sy 

Thus, by setting the twist equal to zero in equation (l) and sub- 
stituting from equation (A50) for Mxy^ elastic constant Cj^ is 
obtained as 

Cjj. = -ctg (A38) 

When the twist is set equal to zero in equation (2) and 

Sx 8y 

as may be determined from equation (A22) is substituted for of 

equation (2), the elastic constant is obtained as 

= Gqtg (A39) 

From the relationships of equations ( 5 ) the remaining elastic constants 
Dxy> <^xy^ become 

Djjy = GJtg5 (A40) 


= (a4i) 

2c^q + 


T = 


Gjtg^ 


(a42) 


Determination of a' and p' in Terms of j and q 

Reference 4 derived expressions for the elastic constants of 
integrally stiffened plates in which the effectiveness of the ribs for 
resisting twisting and shearing is expressed in terms of two param- 
eters a' and P'. The coefficient p' represents the part of the 
rib which is effective in resisting shear when this part is considered 
to be flattened out over the length of the element, thus increasing 
the effective flat-plate thickness of the element. The relation between 
P* and the applied shearing force ^xy is then 
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^xy 



(a45) 


Equating as obtained from equation (A22) to that as obtained in 

eqixation (a 45) results in the following relationship between p' and q: 


P' 


q - 1 

Aw/laS'tS 


(A44) 


The coefficient a' represents the height above the midplane of the skin 
at which the centroid of the distributed fractional area of the rib is 
located to produce the required twisting stiffness. By substituting 
from equation (a 43) for P' and from equation (9) for D^y in equa- 
tion (95) of reference 4, a* is determined in terms of j and q: 


a' 



(a45) 
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TABLE I.- VALUES OF q, c, AND 3 ' AS DETEE^NED BY ANALOG C(»{PUTER 


'•w 

*s 

1 

i 

'S 


q 

c 

P' 

ty,/ts = 0 

1 

2 

1 

1.167 

0.084 

0.778 

2 

h 

2 

1.267 

.145 

.623 

3 

6 

3 

1.373 

.219 

.579 

6 

12 

6 

1.695 

.496 

.540 

8 

16 

8 

1.840 

.645 

.489 

12 

24 

12 

2.090 

.955 

.424 

16 

32 

16 

2.370 

1.224 

.399 

ty/ts = 1 



7 

1 

1.023 

0.018 

0.126 



T 

2 

1.023 

.022 

.080 



T 

3 

1.023 

.024 

.054 



5 

1 

1.033 

.025 

.165 



5 

2 

1.033 

.026 

.083 

0 


5 

3 

1.033 

.026 

.055 



3 

1 

1.058 

.043 

.174 



3 

2 

1.058 

.044 

.087 



3 

3 

1.058 

.045 

.058 



1 

1 

1.240 

.152 

.240 



1 

2 

1.241 

.156 

.121 



1 

3 

1.242 

.158 

.081 



9 

1 

1.070 

.051 

.441 



9 

2 

1.070 

.053 

.259 



9 

3 

1.071 

.053 

.186 



7 

1 

1.092 

.064 

.451 



7 

2 

1.094 

.070 

.274 



7 

3 

1.094 

.070 

.192 

1 


5 

1 

1.138 

. 0 ^ 

.486 



5 

2 

1.138 

.097 

.284 



5 

3 

1.138 

.098 

.201 



3 

1 

1.250 

.143 

.525 



3 

2 

1.250 

.158 

.309 



3 

3 

1.250 

.160 

.219 



7 

2 

1.257 

.165 

.616 



7 

3 

1.257 

.174 

.459 



7 

4 

1.257 

.174 

.367 

2 


‘ 5 

2 

1.401 

.232 

.687 



5 

3 

1.401 

.245 

.512 



5 

4 

1.401 

.245 

.408 



' 9 

3 

1.343 

.245 

.451 



9 

4 

1.343 

.246 

.393 



9 

5 

1.343 

.246 

.348 

5 


7 

3 

1.489 

.516 

.^99 



7 

4 

1.489 

.516 

.435 



7 

5 

1.489 

.516 

.386 

6 

13 

6 

1.748 

.565 

.»^53 

8 

IT 

8 

1.910 

.716 

.437 

12 

23 

12 

2.155 

1.077 

.392 

16 


16 

2.420 

1.361 

.573 



'>s 


q 

c 

P' 

V*S ‘ 2 



8 

1 

1.076 

0.065 

0.304 



8 

2 

1.077 

.075 

.154 



8 

3 

1.078 

.079 

.104 



6 

1 

1.105 

.085 

.515 



6 

2 

1.106 

.101 

.159 



6 

3 

1.108 

.104 

.108 

0 


4 

1 

1.166 

.129 

.352 



4 

2 

1.173 

.153 

.173 



4 

3 

1.173 

.156 

.115 



2 

1 

1.463 

.285 

.465 



2 

2 

1.500 

.340 

.250 



2 

3 

1.500 

.542 

.166 



10 

1 

1.119 

.093 

.490 



10 

2 

1.123 

.126 

.278 



10 

3 

1.123 

.131 

.189 



8 

1 

1.157 

.118 

.517 



8 

2 

1.161 

.158 

.291 



8 

3 

1.161 

.164 

.200 

1 


6 

1 

1.227 

.159 

.561 



6 

2 

1.231 

.213 

.313 



6 

3 

1.233 

.220 

.217 



4 

1 

1.409 

.246 

.673 



4 

2 

1.433 

.331 

.391 



4 

3 

1.440 

. 3^1 

.270 



10 

2 

1.245 

.161 

.462 



10 

3 

1.245 

.164 

.535 



10 

4 

1.245 

.165 

.263 



8 

2 

1.327 

.205 

.492 

2 


8 

3 

1.327 

.204 

.558 



8 

4 

1.327 

.207 

.281 



6 

2 

1.489 

.270 

.551 



6 

3 

1.489 

.273 

.401 



6 

k 

4 

1.489 

.275 

.515 



12 

4 

1.502 

.413 

.405 



12 

5 

1.502 

.413 

.558 




6 

1.502 

.417 

.519 

4 


10 

4 

1.670 

.503 

.451 



10 

5 

1.670 

.504 

.397 



10 

6 

1.670 

.506 

.355 



14 

6 

1.850 

.654 

.423 

6 


■14 

7 

1.850 

.661 

.595 



14 

8 

1.850 

.663 

.569 

8 

18 

8 

1.957 

.840 

.396 

12 

26 

12 

2.185 

1.180 

.559 

16 

34 

16 

2.455 

1.545 

.349 

. 
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TABLE I.- VALUES OF q, c, AND AS DETERMINED BY ANALOG COMPUTER - Concluded 


fH 

5 a 

■^s 

qi 

c 

3 ’ 

V*s = 



6 

1 

1.343 

0.250 

0.515 



6 

2 

1.400 

.380 

.300 



6 

5 

1.424 

.456 

.212 

0 


6 

k 

1.426 

.485 

.160 



k 

1 

1.671 

.395 

.671 



k 

2 

1.852 

.606 

.426 



k 

3 

1.911 

.718 

.304 



4 

4 

1.917 

.768 

.229 



10 

2 

1.509 

.397 

.524 



10 

3 

1.526 

.461 

.383 



10 

k 

1.529 

.489 

.299 

2 


10 

5 

1.551 

.506 

.244 




8 

2 

1.729 

.497 

.600 




8 

5 

1.757 

.576 

.441 




8 

4 

1.765 

.613 

.345 




8 

5 

1.770 

.635 

.284 



12 

4 

1.815 

.682 

.427 

4 


12 

5 

1.816 

.698 

.365 



12 

6 

1.818 

.713 

.318 



12 

7 

1.820 

.713 

.282 



16 

6 

1.950 

.855 

.386 

6 


16 

7 

1.950 

.860 

.350 



16 

8 

1.950 

.865 

.320 



16 

9 

1.950 

.865 

.296 



20 

8 

2.070 

.995 

.356 

8 


20 

9 

2.070 

1.002 

.515 


20 

10 

2.070 

1.006 

.296 



20 

11 

2.070 

1.007 

.280 

12 

28 

12 

2.290 

1.420 

.330 

16 

56 

16 

2.540 

1.880 

.317 


Ik 

5 a 

5 k 


c 

3' 




9 

1 

1.126 

0.107 

0.379 



9 

2 

l.l 4 l 

.164 

.212 



9 

3 

1.143 

.169 

.143 



9 

4 

1.143 

.171 

.107 



7 

1 

1.172 

.141 

.401 



7 

2 

1.190 

.212 

.222 



7 

3 

1.192 

.220 

.150 



7 

4 

1.192 

.222 

.112 



5 

1 

1.263 

.202 

.438 



5 

2 

1.289 

.299 

.241 



5 

5 

1.295 

.510 

.164 



5 

4 

1.295 

.521 

.123 



3 

1 

1.6o4 

.361 

.604 



3 

2 

1.676 

.518 

.358 



3 

5 

1.714 

.550 

.238 



5 

4 

1.714 

.558 

.178 



f 9 

1 

1.228 

.166 

.599 



9 

2 

1.248 

.227 

.347 



9 

5 

1.257 

.251 

.245 



9 

4 

1.259 

.260 

.188 



7 

1 

1.328 

.218 

.670 



7 

2 

1.549 

.293 

.380 

X 


7 

3 

1.359 

.324 

.266 



7 

4 

1.360 

.537 

.203 



5 

1 

1.553 

.331 

.806 



5 

2 

1.582 

.413 

.453 



5 

3 

1.587 

.454 

.511 



5 

4 

1.590 

.471 

.237 



^11 

2 

1.502 

.261 

.430 



11 

3 

1.310 

.289 

.318 



11 

4 

1.311 

.299 

.249 



11 

5 

1.311 

.301 

.205 



9 

2 

1.395 

.319 

.460 

2 


9 

3 

1.406 

.352 

.341 



9 

4 

1.408 

.366 

.268 



9 

5 

1 . 4 o 8 

.369 

.219 



7 

2 

1.573 

.409 

.519 



7 

3 

1.591 

.452 

.386 



7 

4 

1.595 

.470 

.303 



7 

5 

1.595 

.475 

.249 



'13 

4 

1.551 

.490 

.380 



13 

5 

1.557 

.499 

.331 



15 

6 

1.562 

.505 

.294 

4 


13 

7 

1.564 

.507 

.263 



11 

4 

1.730 

.580 

.426 



11 

5 

1.732 

.591 

.368 



11 

6 

1.741 

.596 

.328 



11 

7 

1.743 

.597 

.294 

c 

j 

fl 5 

6 

1.899 

.740 

.403 

D 

1 

[15 

7 

2.020 

.750 

.370 

8 

fl 9 

8 

2.020 

.930 

.385 

|l 9 

9 

2.020 

.940 

.377 

12 

27 

12 

2.250 

1.290 

.356 

16 

55 

16 

2.490 

1.710 

.330 
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Figure 1.- Integrally stiffened plate considered. 


Z 



Figure 2.- Cross section of repeating element and boundary conditions on 

stress function (f>. 
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y 


Figure 


Figure 4 . - 




Repeating element with resultant forces necessary to maintain 
imposed deformation. 
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Figure 5*- Analog Field Plotter. 
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(a) Cross section of repeating element. 
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(b) Values of 

Figure 6.- Values of for repeating element wi 

and. = "by . 


wnen 
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bs/^s 

q 




Figure 7 . - Values of 


^s/'ts 


- 2 




'W 


when both — 

ts 


and 


q ^ 

than those values presented in table I. 


ts 


are larg^ 
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b b b 

Figure 8.- Values of c when both — and ^ are larger than those 

ts ts t3 

values presented In table I. 
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